Photoemission spectra of massless Dirac fermions on the verge of exciton condensation 
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Angle-resolved photoemission spectroscopy (ARPES) is a powerful probe of electron correlations 
in two-dimensional layered materials. In this Letter we demonstrate that ARPES can be used to 
probe the onset of exciton condensation in spatially-separated systems of electrons and holes created 
by gating techniques in either double-layer graphene or topological-insulator thin films. 



Introduction — Electron-hole pairs can undergo Bose- 
Einstein condensation forming a macroscopically coher- 
ent neutral fluid, the exciton condensate (EC) 1 4 . EC 
states have been predicted theoretically in a variety of 
different systems. Spatially separated electrons and holes 
located on two nearby layers (like those that can be cre- 
ated in GaAs/AlGaAs semiconductor double quantum 
wells) were considered first some time ago [S]. These 
double-layer exciton condensates (DLECs) possess spon- 
taneous inter-layer coherence and support counter-flow 
supercurrents. The DLEC has been discussed more re- 
cently in the context of quantum Hall bilayers at total 
filling factor ut = 1 [BHS]- The first clear experimental 
evidence for exciton condensation was obtained a decade 
ago in quantum Hall bilayer DLECs [10], long after the 
state was predicted. 

There is an ongoing effort to find new systems which 
display exciton condensation. In addition to the fun- 
damental interest in understanding when this type of 
order occurs in nature, there is a practical interest in 
finding ECs that occur under less extreme physical con- 
ditions, in particular in the absence of strong mag- 
netic fields and, possibly, at non-cryogenic temperatures. 
Very recently exciton condensation has been predicted 
to occur in double-layer graphene [T1TH4] and in three- 
dimensional (3D) gated topological insulators (TIs) [T5] , 
In the TI case spontaneous coherence is established be- 
tween topologically protected single-particle states on top 
and bottom surfaces. Both double-layer graphene and 
TI thin-film systems are described at low energies by two 
massless Dirac fermion (MDF) Hamiltonians which in- 
teract. Since MDF systems are gapless and truly two- 
dimensional (2D), the field-driven carrier densities that 
can be achieved are much larger than in the semiconduc- 
tor case. Weaker dielectric screening and linearly disper- 
sive conduction and valence bands help to increase both 
interaction and disorder energy scales. Finally, MDF 
bands are perfectly particle-hole symmetric, guarantee- 
ing the almost complete nesting between conduction- and 
valence-band Fermi surfaces, which favors the coherent 
state. 
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FIG. 1: (Color online) The spectral function of a massless- 
Dirac-fermion double-layer system on the verge of exciton 
condensation. The spectral function A(k, ui) (in units of 1/£f, 
where ef is the Fermi energy) evaluated at \k\ = fep is plotted 
as a function of ui/ef for different values of T c /T. The single 
peak at the Fermi energy in the high-temperature limit splits 
into two peaks as T -> T c from above. All curves have been 
evaluated by setting T = 0.1 Tf- 



Graphene double-layer systems can be realized by sep- 
arating two graphene layers by a dielectric [TS] (such as 
AI2O3) or by a few layers of a one-atom-thick insulator 
such as h-BN jT?H2"0"] . In both cases inter-layer hybridiza- 
tion is negligible and the nearby graphene layers are, 
from the point of view of single-particle physics, isolated. 
Topologically-protected massless Dirac 2D electron sys- 
tems (MD2DESs) appear automatically [21] at the top 
and bottom surfaces of a 3D TI thin film. 3D TIs in 
a slab geometry offer two surface states that can be far 
enough apart to make single-electron tunneling negligi- 
ble, but close enough for Coulomb interactions between 
surfaces to be important. 

We report on an analysis undertaken in support of the 
experimental quest for ECs with MDF character. Be- 
cause of the potential value of probes which can detect 
incipient order while still above the critical temperature 
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T c , we examine the influence of pairing fluctuations on 
the ARPES spectra that are commonly used to probe 
two-dimensional systems like graphene. Fluctuation cor- 
rections to thermodynamic and transport observables 
have previously been extensively studied in superconduc- 
tors [52]. General considerations based on the uncer- 
tainty principle allow one to conclude that, like Cooper 
pairs in superconductors, excitons have a finite lifetime 
above T c which is inversely proportional to the distance 
from the critical temperature. Exciton fluctuations will 
therefore influence the behavior of the double layer pro- 
viding distinct signatures in the vicinity of the EC insta- 
bility. Mink et al. |23j . for example, have investigated the 
role of fluctuations in a Coulomb-drag transport setup, 
demonstrating that the drag resistivity grows logarithmi- 
cally upon lowering the temperature towards T c . 

Experimental evidence for ECs in double-layer sys- 
tems has so far been largely based on transport mea- 
surements. In the case of quantum Hall systems the 
transport phenomenology [5 HTUIl2~i] is spectacular in sep- 
arately contacted coherent double layers. In this Letter 
we study the effect of fluctuations on a contactless ob- 
servable: the one-body spectral function, which can be 
accessed by ARPES. ARPES, a photon- in electron-out 
experiment [25] . has been established as an extremely 
important tool in studying the quasiparticle dynamics in 
graphene [2B] and provided the first direct evidence for 
conducting surface states in 3D TIs [27] . 

The quantity measured in ARPES is the one-body 
spectral function A(k,uj). A normal Fermi liquid is 
characterized by a spectral function that, for momenta 
|fc| close to the Fermi wave number fcp, has a single 
Lorentzian peak as a function of frequency with a width 
that vanishes as (k — kp) 2 [23]. Below we demonstrate 
that the spectral function and its integral over k, i.e. 
the density-of-states (DOS) u(ui), show precursors of the 
DLEC phase of unhybridized MD2DESs at temperatures 
larger than the transition temperature. The peak in 
A(k, w)|i fe i_ fcF broadens as one approaches the transition 
and splits into two peaks, as can be seen in Fig. [T] At 
the same time v{uS) develops a strong dip at the Fermi 
energy in contrast to the featureless behavior displayed 
by the non-interacting system (Fig. [3]) . 

MDF Hamiltonian and exciton condensation — We con- 
sider a double layer hosting MDFs whose top (bottom) 
layer is p-doped (n-doped). The two layers are coupled 
only via the Coulomb interaction between the carriers 
and electron tunneling is completely suppressed. Fol- 
lowing Rcfs. [TTI - H5] . we include in our calculations only 
the conduction band (upper Dirac cone) for the n-doped 
layer and the valence band (lower Dirac cone) for the p- 
doped layer. This is justified since the temperatures we 
consider are sufficiently small and the screening lengths 
sufficiently large, so that the energetically remote bands 
can, in first instance, be considered as inert. The band 
dispersions can be approximated by £t>(fc) = u|fe|— £f and 
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FIG. 2: a) The diagram considered for the self energy: the 
upper line refers to the top (hole) layer while the lower line 
to the bottom (electron) layer, b) The ladder approximation 
for the effective interaction V e ff - Eq. |2|. The thin wavy line 
is the screened interlayer interaction Vo{k\,k,2,K)\ K is the 
center-of-mass momentum and p is an integration variable. 



£t(fc) = — v\k\ + ep = — £,b(k) for the bottom (b) and top 
(t) layers, respectively. Here v is the Dirac velocity, which 
is 10 6 m/s for double-layer graphene and ~ 5 x 10 5 m/s 
for the surface states of typical known TI thin films [29] . 
We set h = 1 throughout this work and include only 
inter-layer interactions, assuming that intra-layer inter- 
actions can be taken into account by renormalizing the 
bare parameters of the two separate-layer Hamiltonians. 

We therefore consider the following effective Hamilto- 
nian: 

fc,i— b,t 

+ v o(ki,k 2 , K)cl M _ K cl k2 c hM c tM -K , 

(1) 

where c^fc {c\ k ) is the annihilation (creation) operator 
for an electron in the z-th layer with momentum k. The 
bare (repulsive) interaction Vo(ki, k 2l K) depends on the 
three wave vectors fci, k 2 , and K (and not only on 
l^i — k 2 \, which is the momentum transfer) because of 
the chirality factors [30], which encode the suppressed 
backscattering property of MDFs. After a particle-hole 
transformation c t .fc — > hi k in the top layer one imme- 
diately sees that the second term in Eq. describes 
pairing between an electron with momentum fci in the 
bottom layer and a hole with momentum — ki +K in the 
top layer, interacting via an attractive potential equal to 

-Vo(fci,*a,-K0- 

The EC state is signaled by a pole in the finite- 
temperature many-particle scattering amplitude V g 
which appears at a critical temperature T = T c when 
the center-of-mass momentum K of a pair and the "con- 
jugate" energy Q — £t{K — fci) + £b(fci) vanish. In 
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Fig. [2] we show the diagrammatic representation of the 
electron self energy [panel a)] and the many-particle 
scattering amplitude in the ladder approximation [panel 
b)]. The ladder sum in panel b) is analytically in- 
tractable for a generic two-body potential Vb(fci, k 2 , K). 
To make progress, we approximate [53] Vb(fci, k 2 , K) by 
a momentum-independent constant U. (As discussed in 
Ref. [23] , close to T c the chirality factors do not play any 
qualitatively significant role.) This approximation is rea- 
sonable since it relies on screening, which is robust when 
the electron and hole densities in the two layers are large. 
Substantial numerical work is needed to transcend this 
simple approximation and, in any event, a brute-force 
numerical approach would mask the elementary under- 
standing of the physical effect we want to explain here. 

When the bare interaction does not carry momentum 
labels, the ladder sum in Fig.|2]is a geometric series which 
can be summed to yield 



V eff {K,n u ) = 



U 



l-UE(K,Sl v ) ' 



(2) 



Here f2„ = 2ixv/[3 are bosonic Matsubara frequencies 
with /3 = ObIT 1 , 
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is a pairing susceptibility, A is the 2D electron system 
area ,and np(E) — l/[l+exp(f3E)] is the Fermi-Dirac dis- 
tribution. The pairing susceptibility evaluated at K = 
and = diverges logarithmically at zero temperature, 
and therefore the denominator of Eq. ^ must vanish at 
some non-zero temperature T c no matter how small the 
electron-hole attraction. This is the mean-field pairing 
transition temperature |31j . As temperature decreases 
towards T c , pairing fluctuations first lead to an enhance- 



ment of V c g for small K and VL and ultimately at T = T c 
to a pole at fl = and K = 0. 

Exciton fluctuations and ARPES spectra. - The 
anomalous properties of the effective interaction V e s on 
approaching T c lead to observable consequences [22j . 
Here we concentrate on the spectral function of the bot- 
tom (n-doped) layer Ab(k,ui) = — 7r _1 Sm G^(k,u>), 
which can be calculated from the retarded self-energy 
S^(fc,oj). As is customary [22], we proceed by first eval- 
uating the self-energy on the imaginary-frequency axis 
and then analytically continue it to real frequency. Ac- 
cording to the diagram shown in Fig. [2^), the self-energy 
is given by 



E b (k,uj n ) 



1 

~JA 
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V eS (K, S}„)G[ 0) (k - K, uj n - fi„) 



(4) 

where uj n = 2ir{n + 1/2) //3 are fermionic Matsubara fre- 
quencies and g[° (fc,w„) = [iuj n — Ct(^)] _1 is tne free- 
electron Green's function in the top layer. In order to 
calculate the self-energy from Eq. Q we approximate 
V e s by a form valid [22] at small K and f2 where fluctu- 
ations are strongest: 



-iO 



8/(7^0) 



(5) 



where a(T) = 8 log(T/T c )/(irf3) measures the distance 
from the critical point, £ = [7C(3)w 2 ^ 2 ]/(327r 2 ) is the co- 
herence length, and vq — £-f/(2itv 2 ) is the single-particle 
DOS at the MDF Fermi energy. Using Eq. ^ for the 
retarded fluctuation propagator we can then perform the 
sum over the Matsubara frequencies in Eq. Q. After 
analytic continuation we obtain the following expression 
for the imaginary part of the retarded self energy: 



3mE^(Jb,w) 



d 2 q 



u + 6b(<?) 



tt/^o J (2^) 2 ( w + Cb(g)) 2 + 7(fc-q) 5 



tanh 



coth 



(6) 



where 1 (q) = a(T) + 8£, 2 q 2 /(tt(3). 

Once the imaginary part of the self-energy is known, 
the real part can be calculated via the Kramers-Kronig 
relation CM |35] 



^e^(k,uj) = 3te^(k,oo)+V J * 



dv 3toE^ (fe, v) 



(7) 



Since 3?eEj^ (fc, 00) depends weakly on momentum k near 
kp and on temperature it can be absorbed by measur- 
ing energies relative to the Fermi energy. In practice 
we compute the self energy by numerical integration of 



Eq. ^ and Eq. Q, neglecting 5fte£^(fc,oo) in the lat- 
ter. In the weak-coupling limit we find Ab(k,uj) to be 
a delta function peaked at some frequency Aui. Finally, 
we calculate Ab(k, uj) for various values of T/T c by using 
5fteS^(fe, 00) = —Aui. The parameter we change in order 
to approach the phase transition is thus the interaction 
strength U or, equivalenty, T c . 

Numerical results — In Fig. [T] we plot the sum of the 
spectral functions of top and bottom layers, A(k, uj) — 
Ab(k, uj)+At(k, uj). Due to the perfect particle-hole sym- 
metry of our model, the spectral function of the top (hole- 



4 



3.0 
2.5 
2.0 
3. 1.5 
1.0 
0.5 
0.0 



..... t/t c = 2.oo t\ :a 






T/T c = 1.33 1 V..''W 

\ ^ / 




---- T/T c = 1.11 \ / 




- — t/t c = i.oi vy 




• • • • tunneling 

i i i 


i 



-1.5 



1.0 



-0.5 



0.0 

ui/e F 



0.5 



1.0 



1.5 



FIG. 3: (color online) The density-of-states v [ui) (in units of 
the non-interacting value vq) as a function of ui (in units of 
and measured from the Fermi energy ep). The color coding 
and labeling is identical to Fig. [I] The filled (yellow) circles 
denote the convolution on the r.h.s. of Eq. (J9| evaluated at 
T/T c = 1.01. 



doped) layer At{k,uj) is given by At(k,u>) — Ab(k,—uj). 
In Fig. [I] we illustrate the dependence of A(k, ui) on fre- 
quency ui for |fe| = hp. This plot refers to a fixed tem- 
perature, T = 0.1 Tp, and to different values of T/T c . 
We note that the single Lorenzian peak of the spectral 
function becomes broader and broader as one approaches 
the phase transition and eventually splits into two sep- 
arate peaks very close to the phase transition. Note 
that the splitting between the two peaks is very large, 
~ 5&bT c at T/T c = 1.01. The two-peak structure resem- 
bles the behavior of the spectral function expected below 
the transition temperature. From Bogoliubov theory one 
indeed expects two delta-function peaks at frequencies 
fl± = ±|A| where A is the excitation spectrum gap. 

From the spectral function one can also calculate the 
DOS v{u): 



v{u>) 



d 2 k 



A{k,uj) 



(8) 



We performed the integral numerically after evaluating 
the integrand for momenta up to |fe| = "ik-p. Our results 
are shown in Fig. [3j where, for the sake of comparison, 
we have also plotted the non-interacting result. We note 
that as one approaches the phase transition more and 
more spectral weight is shifted from the Fermi energy 
(a; = 0) to higher energies. The resulting dip in the 
DOS around ui = can be seen as a precursor of the 
gap in the excitation spectrum that opens up only in the 
broken-symmetry phase. 



In recent decades the DOS of electronic systems has 
been measured by using scanning tunneling microscopy 
(STM), which probes locally to mitigate the influence of 
inhomogeneities and some types of disorder |36| . For ex- 
ample the zero bias anomaly of disordered metals [37J has 
been successfully explored in this way. The differential 
tunneling conductivity is normally well approximated by 



(v\ dI 





dn F (E)~ 


I — oo 


BE 



v{E + eV) dE . (9) 



In Eq. ([9| / and V are respectively the tunneling cur- 
rent and the bias voltage between the STM tip and the 
sample. When v (E) is smooth on the scale of k^T the 
distribution function on the r.h.s. of Eq. ([£]) acts like 
a delta function so that otunOO oc v(V) [37]. The sit- 
uation is very different in the case of STM studies of 
fluctuating superconductors above T c . In this case the 
DOS can have structure that is sharp on the scale of the 
temperature. The tunneling conductivity <7tun(^0 then 
has a broad pseudogap structure (with a maximum at 
eVmax = wksTc), resembling the BCS gap below T c [35] , 
This occurs in spite of the predicted [39 , 40 j sharp singu- 
larity in the fluctuation contribution to the DOS close 
to T c . The tunneling conductivity is much smoother 
as a function of temperature or voltage (<Tt U n(0)/crN oc 
In [T C /(T — T c )] where <tn is normal conductivity of the 
tunnel junction) than DOS is as a function of energy. 
This distinction ultimately reflects the property that the 
energetically averaged DOS is not altered by order [4"T] . 
Similar considerations apply to the EC fluctuation con- 
tribution to the tunneling conductivity as illustrated in 
Fig. [3] (circles) where we have plotted the r.h.s. of Eq. ^ 
at T/T c = 1.01. 

In summary, we have demonstrated that ARPES and 
STM can be used to probe the onset of exciton condensa- 
tion in spatially-separated systems of electrons and holes 
created by gating techniques in double-layer graphene 
and topological-insulator thin films. The advantage of 
these experimental probes is that one does not need to 
realize independent contacts to two closely-spaced Dirac 
layers. These probes can be particularly useful in the 
context of topological-insulator thin films where separate 
contacting is precluded by hybridization between top and 
bottom at the sample edges. 
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